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RESEARCH MEMORANDUM 

THE DESIGN OF VARIABLE MACH NUMBER ASYMMETRIC SUPERSONIC NOZZLES 
BY TWO EROCEDURES EMPLOYING INCLINED AND CURVED SONIC LINES 
By Clarence A. Syvertson and Raymond C. Savin 

SUMMARY 

Two theoretical procedures are developed for designing asymmetric 
supersonic nozzles for which the calculated exit flow is nearly uniform 
over a range of Mach numbers. One procedure is applicable at Mach num- 
bers less than approximately 3- This approach yields, without itera- 
tion, a nozzle for which the calculated exit flow is uniform at two 
Mach numbers and, with proper design, is nearly uniform at Mach numbers 
between, slightly above, and slightly below these two. The use of an 
inclined and curved sonic line is an essential feature of this approach. 
The second procedure requires iteration and is used for designs at Mach 
numbers exceeding 3* Although it is not a necessary feature, an 
inclined and curved sonic line is also used in this procedure. In both 
approaches the flow field downstream of the sonic line is determined 
using the method of characteristics. 


INTRODUCTION 

The asymmetric, sliding-block— type nozzle was first shown by Allen 
(reference l) to be an effective means for producing reasonably uniform 
supersonic flow over a range of Mach numbers. Allen indicated that the 
principal advantage of this type of nozzle was its structural and 
mechanical simplicity; however, he noted that the asymmetric nozzles 
which were investigated had shown undesirable horizontal and vertical 
static— pressure gradients at some operating speeds. This observation 
was later confirmed by the results of Frick and Olson (reference 2). 

An investigation was therefore undertaken to develop systematic proce- 
dures for designing asymmetric nozzles which would produce exit flows 
of the desired uniformity for wind-tunnel testing over a given range of 
Mach numbers. The purpose of the present paper is to describe the two 
most suitable design procedures developed to date at the Ames Labora- 
tory. Many features of these procedures have resulted primarily from 
experience gained in the design of 18 nozzles, 8 by the first procedure 
to be presented, and 10 by the second procedure. 
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SYMBOIS 


local velocity of sound 
height 

Mach number, ratio of local velocity to local velocity of 
sound 

static pressure 
dynamic pressure 

radius of curvature of streamline 

dimensionless perturbation velocity component parallel to 
x axis l — 1 


dimensionless perturbation velocity component parallel to 
v 


y axis l — 
J V a* 


velocity component parallel to x axis 
velocity component parallel to y axis 


dimensionless resultant velocity 


—2 __ 2 

U + V 


a* 


rectangular coordinates 


Mach angle (^arc sin 


i) 

M / 


ratio of specific heats (1.400 for air) 


stre am static-pressure coefficient 


(¥>) 


' "*0 ' 

angle between direction of flow and some reference axis 


Subscripts 


lower wall 
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M conditions at a specified Mach number 

N measured normal to a surface 

o initial or reference conditions 

s conditions along a streamline 

U upper wall 

Superscripts 

* conditions where the Mach number is unity 


ANALYSIS 

Two Basic Features of the Variable Mach Number Asymmetric Nozzle 


Before discussing specific procedures for designing two-dimensional 
asymmetric nozzles, it is useful to consider from a q ual itative point 
of view the reasons why such nozzles are capable of producing uniform 
exit flows at more than one Mach number. It is generally accepted 
that if a supersonic nozzle is to produce uniform exit flow at more 
than one Mach number, the shape of the nozzle passage must be different 
for each exit Mach number. For this reason, most variable Mach number 
wind tunnels have employed either a series of interchangeable nozzle 
walls (blocks), each of fixed geometry, or flexible nozzle walls. With 
both of these solutions the contours of the nozzle walls are changed to 
give the changes in passage shape required for variation of the exit 
Mach number. 

In contrast to these solutions, Allen’s proposal is that the con- 
toured nozzle walls be rigid and that they be moved relative to each 
other so as to form asymmetric passages of various shapes for producing 
flows at different exit Mach numbers . An asymmetric nozzle of this 
type is shown in figure 1. The solid lines show the passage shape for 
a high Mach number exit flow; the dashed line shows the lower rigid 
wall moved to a new position to produce a passage shape for a low Mach 
number exit flow. The lower wall is moved in a direction parallel to 
the walls of the test section and the ratio of minimum section area to 
test section area is thus changed while the test section area is 
maintained constant. This first feature of the sliding-block asymmetric 
nozzle makes it a practical device for varying the exit Mach number. 



k 


NACA RM A51A19 


It is also clear from figure 1 that certain portions of the walls 
are included in the supersonic flow region at the high Mach number 
setting that were included in the less critical subsonic flow region at 
the low Mach number setting. They may, therefore, be shaped to make the 
remaining downstream portions of the walls (designed to produce uniform 
exit flow at the lower Mach number) consistent with uniform exit flow at 
the higher Mach number. This second feature of the asymmetric nozzle 
makes it possible to produce uniform exit flows at more than one exit 
Mach number. As is discussed later, it is this feature which allows the 
development of an explicit procedure for designing asymmetric nozzles 
operating at Mach numbers less than approximately 3 • 


Determination of Flow Characteristics in 
Asymmetric Nozzles 


The problem of determining the characteristics of the flow between 
the walls of an asymmetric nozzle can be divided into three parts. 

Figure 2, a schematic diagram of the construction of the flow field in 
a typical nozzle, will be used to illustrate the different parts of the 
problem. Part 1 is the determination of the flow conditions in the 
region of the minimum section where supersonic flow first occurs. 

Line AB in figure 2 represents the location of the points where the 
initial flow conditions are calculated. Proceeding downstream, part 2 
is the determination of the flow characteristics from these initial 
points to the region just upstream of the nozzle exit (area ABDC in 
fig. 2). Part 3 is the construction of the flow field which must exist 
in the region just upstream of the nozzle exit in order to have uniform 
exit flow. The region of the nozzle affected by this part of the 
problem is represented by the area CDE. 

It will be seen later in the analysis that the procedure developed 
herein for Mach numbers less them 3 requires the nozzle walls to be 
curved in the region of the minimum section. This curvature, of the 
walls will cause the sonic line to be inclined with respect to the flow 
direction and generally curved. Part 1 of the problem of determining 
the flow characteristics in an asymmetric nozzle is the determination 
of the flow in the region of the inclined and curved sonic line. The 
solution of this problem is discussed in appendix A. 

Throughout this study the method of characteristics will be 
employed to construct the region of the supersonic flow field associated 
with part 2 of the problem of determining the flow in asymmetric 
nozzles. The characteristic solution of two-dimensional supersonic 
flows was originally developed by Prandtl and Busemann (reference 3) . 
This solution may be accomplished by numerical or graphical processes, 
or a combination of these two. Of the alternative approaches, the 
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grapMcal method was chosen as most suitable for the nozzle— design 
problems considered here. The reasons for this choice are (l) the net 
size can be easily varied during the construction of the flow field; 
and ( 2 ) the designer is enabled to visualize graphically the velocity 
diagram of the flow field. 


Sauer (references 4 and 5) has developed a graphical method for 
axially symmetric flow fields in which the characteristic lines in the 
holograph plane are constructed in a step— by— step process. An adapta- 
tion of this method to two-dimensional supersonic flows is employed in 
this paper. The hodograph on which the nozzle flow fields are 
diagrammed is constructed using the polar coordinates W and 0 and, 
for convenience, the hodograph is marked with lines of constant and 

0. (A W— scale of 0.02 to the inch is used.) For an irrotational, 
two-dimensional flow field, the characteristic lines in the hodograph 
plane all have the same shape. This property permits the use of a 
characteristic tenqdate which greatly increases the speed of the calcu- 
lations (i.e., the construction of the Mach net). The shape of the 
template can easily be determined since the relationship between Mach 
angle and flow direction along the characteristic lines is the same as 
for the expansion about a corner of a uniform, two-dimens ional , super- 
sonic stream. This relationship, given by the familiar Pr andt 1-Meyer 
equation (see, e.g., reference 6 ), 


0 - 0 * 




tan ( 90 °-a^) - ( 90 °- a^) 


is plotted on the hodograph and the template is made from the curve. 
All the characteristic nets to be discussed are constructed using the 
hodograph and template Just described. 

The requirement for uniform exit flow in a nozzle is satisfied if 
the entire flow field downstream of the nozzle exit can be represented 
by a single point in the hodograph diagram. This point is unique 

since it represents the design flow parameters 0 and a for the 

M 

test section. The effect of this requirement on part 3 of the problem 
of determining the flow characteristics in a nozzle is that area CDE 
in figure 2 will be represented in the hodograph diagram by a single 
second-family characteristic line . 1 This characteristic line passes 
through the unique point mentioned above. The flow in region CDE will 
be of the Prandt 1-Meyer type. A discussion of this type of flow is 
presented in appendix B. 


X A second— family characteristic line is a Mach line negatively inclined 
with respect to the flow direction. 
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Design Procedure for Mach Numbers Less Than 
Approximately 3 


The following procedure for designing variable Mach number asym- 
metric nozzles is the more direct of the two to be presented. This 
procedure consists of three major steps. These steps are as follows: 
First, the contours of the downstream portions of the walls are deter- 
mined to produce uniform exit flow at the lower design Mach number; 
second, the flow field that must exist between these wall portions for 
uniform exit flow at the higher design Mach number setting is deter- 
mined; and third, the portions of the walls upstream of those considered 
in step one are shaped so as to produce the flow field determined in 
step two. These three steps ensure that the nozzle will produce 
theoretically uniform exit flow at the two design Mach numbers. It is 
required, however, that the exit flow be nearly uniform over a range of 
Mach numbers from slightly below the lower to slightly above the higher 
design Mach number. Experience gained in this investigation has indi- 
cated that specifying continuous pressure gradients in the nozzle at 
the design Mach numbers and continuous variation of the second deriva- 
tives of the nozzle wall contours tends to minimize the calculated flow 
irregularities over the Mach number range. These specifications, there- 
fore, are embodied in both design procedures. The uniformity of the 
flow at other than the two design Mach numbers should, of course , be 
investigated. 


It has also been found for the proposed design procedures that in 
order to minimize the flow irregularities at off— design Mach numbers it 
is desirable to have the position of the minimum section move continu- 
ously upstream along both walls as the exit Mach number is increased. 

The continuous movement will be accomplished if the wall contours 
satisfy certain conditions in the region where the mini mum section 
moves. These conditions are: First, the wall contours should have 
single curvature (no inflection points); second, they should have finite 
radii of curvature ; and third, the local radius of curvature of the 
lower wall should always be sufficiently less than the local radius of 
curvature of the upper wall to ensure that in a given position the 
nozzle has only one minimum section. 


The upper and lower design Mach numbers and the test-section 
height are first selected for any design undertaken. It follows from 
the discussion of the uniformity of the exit flow at off— design Mach 
numbers that the two design Mach numbers may be selected somewhat 
within the extremes of the desired range of operation of the nozzle. 

A typical low Mach number construction, representing the first step for 
such a design, is shown in figure 3 . For the first part of this con- 
struction the flow parameters 6 and at a series of initial 

points just downstream of the sonic line are determined from the tran- 
sonic equations given in appendix A. These initial points are shown 
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along line A— B. Use of "the equations in appendix A requires that the 
radii of curvature of both vails in the region of the m-in-iTmim section 
he known. The radius of curvature of the lower wall is made large 
compared to the mi ni .mum section height (of the order of 10 mi n-? -rmnn 
section heights) to permit an accurate application of these equations. 
The upper-wall radius of curvature at the mini mum section is selected 
in accordance with the discussion of appendix B to give Prandtl— Meyer 
type flow in the region .just upstream of the nozzle exit. (See 
equation (B2) . ) To satisfy the conditions of appendix A, both radii 
are kept constant a short distance (of the order of one— half the mini- 
mum section height) upstream and downstream of the min-iTm-i-m section. 

The upper and lower walls are extended upstream to points J and L, 
respectively, to fulfill this requirement. 

For the second part of this construction the characteristic net is 
constructed downstream of the initial points putting in one second- 
family Mach line at a time. Line C— D in figure 3 is a typical second- 
family Mach line. The radius of curvature of the upper wall is held 
constant until the flow along that wall becomes of the Prandtl-Meyer 
type. This occurs at point E of figure 3(a). The curvature of the 
lover vail is held constant until the design Mach number is approached 
on that surface. 3 The curvature is then gradually decreased, becoming 
zero vhere the design Mach number is attained (point G in figure 3( a ))* 4 
When the characteristic net is complete, the shape of the dovnstream 
portions of both vails is determined. The first major step in the 
design procedure is now complete. 


2 The flow along the upper vail becomes of the Prandtl-Meyer type when 
the line representing that vail in the hodograph diagram, line ACE 
in figure 3(b), becomes tangent to a second-family characteristic 
line such as lineE-H. Experience has indicated that this tangency 
will occur if the upper-wall radius of curvature is selected by use 
of equation (B 3 ) and if the other su gg estions are followed. 

3 The curvature of the lower vail is maintained constant over most of 
the length of this construction since this results in the shortest 
nozzle without increasing the curvature. Increasing curvature in the 
streamwise direction has been found difficult to handle in later 
constructions . 

4 It is difficult to establish a definite rule for the length of the 
portion of the lower wall along which the curvature is decreased to 
zero. A length of one-half to one test-section heights has been 
found satisfactory. 
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The flow field that must exist between these portions of the walls 
for uniform exit flow at the higher Mach number setting will now be 
determined. Maintaining the test-section height constant, the upper— 
and lower-wall portions are rotated to make 0=0 in the test section. 
Then holding the upper wall stationary, the lower wall is translated 
horizontally until it is in the higher design Mach number setting. A 
satisfactory criterion for this position is that the downstream ends of 
the contoured portions, points G and H in figure 4(a), be alined along 
the sams first— family Mach line inclined at the Mach angle corresponding 
to the higher design Mach number. Specifying uniform flow downstream 
of this Mach line, the characteristic net can be constructed between the 
portions of the walls in the reverse of the flow direction. Such a 
construction is shown in figure 4. The construction is continued to the 
upstream end of the upper-^wall portion, point J. In general, the 
termination of the construction on the lower wall, point K of figure 4(a), 
will lie downstream of the end of the predetermined portion, point L; if 
not, the lower wall can be extended at constant curvature. This charac- 
teristic net construction determines the flow conditions that must exist 
between the two wall portions if the exit flow is to be uniform at the 
higher design Mach number. 

The next problem is to determine the shape of the portions of both 
walls between the new minimum section and the line J— K of figure 4. 

This shape must be such as to produce the same variation of Mach angle 
and stream angle along the Mach line J— K as was obtained by the 
characteristic net construction outlined in the preceding paragraph. 

The radius of curvature of the lower wall at the minimum section is 
selected in the same manner as for the lower design Mach number. The 
height of the minimum section may be calculated by use of the one- 
dimensional area— ratio relationship since the curvature of the sonic 
line does not significantly influence the mass flow. The radius of 
curvature of the upper wall is selected with the aid of equation (B2) 
to give approximately the same variation from Prandtl-Meyer type flow 
as occurs at the Mach line J— K. 5 This selection need only be an 

5 The difference from Prandtl-Meyer type flow at point K in figure 4 is 
defined here as the difference between d0/da^ along line KG evalu- 
ated at point K and d0/da^ along the second— family character- 
istic line passing through point K. If, for example, | d0/da^j along 

line KG is 10 percent less than along the characteristic line, the 
curvature of the upper wall should be about 10 percent less than that 
calculated for Prandtl-Meyer type flow. The curvature for Prandtl- 
Meyer type flow may be calculated by use of equation (B2). In that 
equation is the radius of curvature of the lower wall at point K 
and hjj is the normal distance between the two walls which may be 
estimated graphically in figure 4(a). 
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estimate. It should "be pointed out that there is no way of know ing in 
advance the location of the new minimum section with respect to the 
test section. For this reason a different coordinate system is used 
for the new construction. 

Following the procedure of the lower Mach number design, the 
initial points are selected along line A*-B l of figure 5* The charac- 
teristic net will be constructed downstream of the initial points until 
the Mach number on the lower wall coincides with the Mach number at 
point K in figure 4( a) . At this point on the lower wall two other 
conditions should also be satisfied. These conditions are: (l) the rate 
of change of Mach angle with flow angle along the lower wall should be 
the same as at point K, and (2) the radius of curvature of the lower 
wall at this point and upstream of the point a distance equal to the 
distance K— L of figure 4(a) should be the same as at point K. These 
conditions may be satisfied by using the follow in g procedure . The net 
is first constructed as far as point P in figure 5, while simultaneously 
the radius of curvature of the lower wall is gradually increased. This 
part of the construction determines line B*— P in the hodograph 
(fig. 5 ( 12 ) )• The point where the Mach number coincides with point K 
will be designated as point M. A line representing the lower wall in 
the hodograph diagram is faired from point P to point M in a manner to 
be tangent to the line B’— P and to give the same d9/do^ along this 

line at point M as the previous construction (fig. 4) had along the line 
K-G at the point K. This curve represents a lower wall which will 
satisfy condition (l) and determines the value of 9 at point M. With 
the value of 6 and the radius of curvature at point M the shape of the 
portion of the lower wall from point P to point M may be chosen to 
satisfy condition (2). Several points are chosen along line P— M and are 
combined with the Mach line Q— P to determine the net in the region QHMP 
in figure 5(a). The shape of the upper wall in this region is the 
streamline through point Q. 

The two constructions in the physical plane (figs. 4(a) and 5(a)) 
are then joined at points K and M and properly oriented by equating 6 
at the two points. This orientation brings the two constructions to 
the same coordinate system. The combined constructions in the physical 
plane are shown in figure 6. The characteristic net is constructed 
between the two portions in the region NKJ and the upper wall contour is 
determined by plotting a streamline connecting the two portions of the 
upper wall. This completes the design of the contoured walls which will 
produce theoretically uniform exit flow at both design Mach numbers. 

The flow in the nozzle should be investigated at several other Mach 
numbers in the design range to determine if any significant variations 
of Mach number or flow inclination occur in the test section. This has 
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not occurred in any of the eight designs attempted. Should any signi- 
ficant variations occur, however, appropriate corrective measures of the 
type discussed in the next section can he taken. 

In theory there is no apparent reason why this design procedure 
should he restricted to two Mach numbers or why it should he limited to 
Mach numbers less than 3- Practical difficulties do limit its applica- 
bility, however. For example, if the higher design Mach number exceeds 
approximately 3, it becomes very difficult to find contours for the 
lower Mach number portions of the walls to which any additional portions 
may he added to give uniform exit flow at this higher Mach number. The 
difficulty develops when the characteristic net is constructed in the 
reverse of the flow direction. Expansion wavelets, which diverge in the 
streamwise direction, converge in this construction. Above M = 3 the 
rate of convergence is such that, in approximately 10 constructions 
attempted, several waves were found to coalesce. When this coalescense 
occurs, no solution at the higher design Mach number is possible using 
the present method. A similar difficulty (note discussion in appendix B) 
develops when the procedure is repeated for a third Mach number below 3* 
These difficulties could possibly be overcome by considerably increasing 
the length of the nozzle. This solution is impractical, however, for 
many applications (e.g., large wind tunnels). 


Design Procedure for Mach Numbers Greater Than Approximately 3 


The procedure for designing variable Mach number asymmetric nozzles 
to operate at Mach numbers greater than 3 employs an averaging process 
and yields nozzles for which the analysis indicates approximately uni- 
form exit flow. The contour is designed first for the highest Mach 
number by a procedure similar to that used for the design of the lower 
Mach number portion of the walls discussed in the preceeding section. 
Figure 7 represents the physical plane of a construction of this type. 

A difference between this design procedure and the one referred to above 
is that the flow along the upper wall does not become of the Prandtl— 
Meyer type before point C. At this point it must, to satisfy the 
requirements for uniform exit flow. In general, it has been found 
desirable to turn the flow between points A and E through an angle no 
greater than the Prandtl-Meyer expansion angle for the lowest design 
Mach number. The reason for this limitation is discussed in the next 
paragraph. 

The contours for the highest Mach number are translated to the 
lowest design Mach number setting and the characteristic net is recon- 
structed. Adjustments are made to the shape of the portion of the upper 
wall between point C and point E of figure 7 to provide uniform exit 
flow at this lowest Mach number. The process is repeated at two 
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intermediate Mach numbers. For the portion C— E, four contours result. 
These results are plotted in the form of slope curves (i.e. , 6 as a 
function of x) and an arithmetic mean of the extremes taken. Exper- 
ience gained in constructing 10 nozzles by this method has shown that 
limiting the angle through which the flow is turned in the manner sug- 
gested in the previous paragraph will minimize the difference between 
the four contours for the portion of the wall between point C and 
point E. 

The mean contour is then substituted for the portion C— E in each of 
the four constructions. The characteristic nets are changed in this 
region and extended into the test section. From these constructions, 
plots are made of the stream angle arid Mach number distributions in the 
test section. Should any unacceptable variations result, they are 
traced back through that particular construction and appropriate adjust- 
ments made to the walls. The effect of these adjustments on the other 
three constructions must then be considered. Through a series of such 
adjustments, nearly uniform flow (stream-angle variations of less than 
±0.25°) can generally be obtained at all four Mach numbers. 

Computation of Ordinates and Boundary-layer Correction 

A first determination of the basic wall contours has now been made; 
however, as was discussed previously, these contours must be examined 
to ensure that the second derivative curves are fair. The slope curves 
of the contours (i.e., 6 as a function of x) are converted to first 
derivative curves (i.e., dy/dx as a function of x) which are differ- 
entiated graphically to obtain the values of the second derivatives. 
Curves are faired through these values and revised first derivatives 
and slopes are obtained from graphical integration of these curves. 
Experience has indicated that the original and revised slopes should 
agree to within ±0.03° at all stations. The process of fairing and 
integrating should be repeated, if necessary, until this criterion is 
satisfied. The reviseu. ordinates are obtained in the usual manner from 
the revised first derivative curves. 

The results of Tucker (reference 7 ) may be used to compute the 
boundary— layer correction for both low and high Mach number nozzles. 

To be consistent with the determination of the uncorrected nozzle con- 
tours, care should be taken to ensure that the boundary— layer correction 
is also fair to the second derivative. 

DISCUSSION 

Many vsxisrfcioiLS of *tlie procedures given for the design of variRhle 
Mach number, asymmetric, supersonic nozzles are possible. The 
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suggestions given result primarily from experience. In regard to the 
fairness of the nozzle walls, something more should be said. As the 
second derivative of the nozzle contour is directly related to the pres- 
sure gradient, continuous variation of the second derivative along the 
nozzle wall tends to ensure continuous variation of the pressure 
gradient. Theoretically, a discontinuity in the second derivative is 
permissible, provided the Mach line (or its reflection) passing through 
this discontinuity also passes through another discontinuity of opposite 
nature and equal magnitude. With asymmetric nozzles of the type under 
consideration, this requirement would not necessarily be satisfied at 
all settings; hence, continuity of the second derivative of both walls 
is specified. It follows, of course, that the curvature of both walls 
goes to zero continuously at the nozzle exit. 

The usefulness of the design procedure for Mach numbers less than 
approximately 3 lies in its direct, noniterative approach and in the 
fact that theoretically 'uniform exit flow in a nozzle can be obtained at 
two Mach numbers. If the flow is uniform at two Mach numbers , it seems 
reasonable to expect the flow to be nearly uniform at intermediate Mach 
numbers and Mach numbers slightly above the higher and slightly below 
the lower design Mach number. Theoretical analyses of the flow at 
several off-design Mach numbers have substantiated this expectation in 
all eight designs completed. 

No experimental results for large-scale nozzles are available with 
which the theoretical predictions of the methods presented in this 
paper may be compared. Some data are available which were obtained in 
a 2- by 2-inch nozzle designed by the procedure for Mach numbers less 
than 3. This nozzle was designed to have uniform exit flow at Mach 
numbers of 1.5 and 2.6. The coordinates for the walls of this nozzle 
are presented in table I. These coordinates do not include a correction 
for boundary layer; however, the coordinates of the nozzle fabricated 
did have a boundary— layer correction.® Static— and pitot— pressure 
surveys were made in the test section and the results of these surveys, 
converted to stream static— pressure coefficient, are shown in figure 8. 
These results indicate that variations in stream static— pressure coeffi- 
cient are small, never exceeding approximately ±0.006 over the Mach 
number range. These variations are of the same order as the errors 
involved in the measurements . The results of the surveys indicate none 
of the large variations in stream static— pressure coefficient that were 
observed by Frick and Olson (reference 2) in the tests which led to the 
present investigation. 

s The coordinates of the walls of this nozzle are not given with correc- 
tion for boundary layer since these coordinates would apply only for 
the Reynolds number at which these tests were conducted. 
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CONCLUSIONS 

Two theoretical procedures are developed for designing asymmetric 
supersonic nozzles for which the calculated exit flow is nearly uniform 
over a range of Mach numbers. One procedure is applicable at Mach 
n umb ers less than approximately 3 and yields, without iteration, a 
nozzle producing theoretically uni form exit flow at two Mach nunibers and 
nearly uniform exit flow over a range of Mach numbers. An essential 
feature of this procedure is the use of an inclined and curved sonic 
line. The other procedure requires iteration and is used for nozzles 
designed to operate at Mach numbers greater than approximately 3- 
Although it is not an essential feature, an inclined and curved sonic 
line is also enqployed in this procedure. 

A 2— by 2— inch nozzle was designed by the procedure for Mach 
n umb ers less than 3 to operate over the Mach number range from 1.5 to 
2.6. Hie results of static— and pitot— pressure surveys in the test 
section of this nozzle indicate a maximum variation of ± 0.006 in the 
stream static— pressure coefficient over the design Mach number range. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. 
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APPENDIX A 


DETERMINATION OF THE FLOW FIELD NEAR THE SONIC LINE 


Sauer (reference 8), among others, has developed equations defining 
the sonic line in flow through the minimum section of a symmetrical 
nozzle. Sauer also developed similar equations defining transonic flow 
about profiles. The equations for the transonic flow about a profile 
may be applied to flow through the minimum section of an asymmetric 
nozzle. 


The origin of the coordinate system is placed at the point where 
sonic velocity occurs on the upper contour, and the x axis is made 
tangent to this contour. (See fig. 9.) By the use of Sauer's develop- 
ment but adding a term in the fourth power of y in the velocity 
potential, the following equations result: 


u = ax - Py + y 2 


a* ♦ 


v = — (3x + xy 


2p 2 + (7+I) a 2 


+ • - • 


<*p + 1+px (^|± a 2 + p*) 


r 7+l 


(Al) 

y 2 + 



(A2) 


(A3) 

(Ah) 


Equations (Al), (A3), and (A4) are the same as Sauer's results. Equa- 
tion (A2) includes one additional term. 


In order to fix the sonic line with respect to the minimum section, 

it is necessary to know x« or x« • (Note fig. 9») It may be 

U JL 

deduced from the geometry of the problem that 




L U 


- x, 


•Mt 


y s 


= (jL. 

\ l+u 


(A5) 


Cl 


x_x Mu 

y=o 
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and 

% R 1 


Setting 


— v ' 

W x=o 

t = y s 

making the necessary substitutions in equations (A5) and (A 6), and 
solving for Xj^ and x^ yields 



(S-u) „ 


y=y s 


\ 


(A 6) 


X M, 


-1 
U 2ay s 


•J (Rl — R lT"ys ) 2 + WsKRuRl + (Rj^- Ry - Ys) 


(A7) 


and 



Equations (Al) through (A8) define the shape of the sonic line, 
provided the local flow velocity is not too different from sonic 
velocity in the region near the sonic line and provided the radii of 
curvature of the lower and upper walls are essentially constant in this 
region. The stream angle and Mach angle at any point near the mini- 
mum section can he calculated by standard methods from the known values 
of u and v. The correction to mass flow due to the curvature and 
inclination of the sonic line was investigated and found to he negli- 
gible for nozzles having moderately curved and inclined sonic lines. 


If the nozzle is symmetrical, then the upper wall is the plane of 
symmetry and it follows that 

3=0 


The previously developed equations for flow in the region of the sonic 
line now reduce to 


7+1 
2 

v = ( 7 +l)a 2 xy + 


/ T-L. p O 

u = ax + -- a y 2 


a/ 


2 — 


iZ+iil a ^y3 
6 

1 1 


7+1 y s R L 
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and 


*m = - ~r ay s 2 

These equations are identical to those developed by Sauer for a 
symmetric nozzle. It should be noted that in this analysis y s is, of 
course, always negative. 
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APPENDIX B 
PRANDTIr-MEYER FLOW 


In the region just upstream of the exit from an asymmetric super- 
sonic nozzle, a special type of flow occurs. This region (ABC in 
fig. 10) is hounded by the curved upper tunnel wall and the Mach line 
passing through the downstream end of this wall. Since the walls are 
straight downstream of the line A— B, any variations in Mach number and 
flow angle along this line will be present in the test region. For 
uniform exit flow the flow conditions along this line, or any other 
first— family Mach line in this general region, must be constant. The 
relationship between Mach number and flow angle along the upper wall 
from point C to point A must therefore be the same as for Prandtl-Meyer 
expansion of a uniform, two-dimensional, supersonic stream about a 
corner. The solution to this flow problem is well known (see, e.g., 
references 9 10) and will not be presented here, but rather the 

application of the solution to the problem of supersonic nozzle design 
will be discussed. 

Consider now a nozzle made up of elements of two Prandtl-Meyer 
stre aml ines as shown in figure 11. It is clear from the figure that 

% - - *m l (bi) 

Assume that the ratio between the normal height from the streamline to 
the effective corner h^ and the radius of curvature of the streamline 
R is a function of the local Mach number only. This will be proven 
subsequently. Replacing b^ and hjj^ with the product of this ratio 

and the proper radius of curvature, equation (Bl) becomes 

% = X (R U - R l) (B2> 

The function 


By combining 


applies only 
streamlines . 
any two Mach 


P is defined as 
P = 


h _ ^m 

r u - r l r 


(B3) 


equations (B2) and (B3) it can be seen that the equation 

P = 0 (B4) 


wnen me rxow is oi me rraaan- meyer type oe tween tne two 
With this type of flow the change in flow angle between 
lines of the first family, such as those shown 
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in figure 11, is the same along either streamline. 

Consider now when 

P >0 (B5) 

If all the quantities in P are treated as fixed with the exception of 
Rj (note equation (B3)) the inequality means that R L is greater than 
for the case of Prandtl-Meyer flow. The change in flow angle along the 
lower streamline between two Mach lines will he less. Some additional 
turning of the stream will he required hy reflected wavelets from the 
upstream portion of the nozzle. In this case the wavelets would he 
expansions . 

Consider the inequality 


P < 0 


(B6) 


By use of the same reasoning as in the previous paragraph, it can he 
shown that reflected compression wavelets are required. 

For a given asymmetric nozzle the ratio hjj/(R[p-RL) will he fixed 
hy the geometry of the nozzle (as determined, for example, at the first 
design Mach number) while h^/R will vary with the Mach number in the 
region of the nozzle exit. Hence, knowing the variation in h^/R it 
will he possible to determine qualitatively the variation in P and 
thus the type and concentration of wavelets that must he reflected in 
the nozzle at various Mach numbers. 


It now remains to determine the manner in which h^/R varies with 
Mach number. The equation of a Prandtl-Meyer streamline in polar coor- 
dinates as given by Sauer (reference 5) is , 


r 


= h. 



(B7) 


The coordinate system is shown in figure 12. The radius of curvature 
of this streamline can he found as a function of h Q and cd. The angle 
cd is related to the local Mach number hy the following equation, also 
taken from reference 5: 


(7+1) - 2 


cos 


M = 



(7-1) 


COS' 


U£ m ) 


(B8) 
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Equation (b 8) can be combined with the expression for the radius of 
curvature in terms of h D and a> to give 


ho (7+1) M 3 

K = 

2 (M 2 — l) 



z±i 

2 


7+1 

2(7-1) 


(B9) 


From one— dimensional area-ratio considerations (reference 6) 


ho 



7+1 

2 

1 + M 2 
2 


7+1 

2(7-1) 


Combination of equations (B9) and (BIO) yields 

= 2 (M g - 1 ) 

B ~ (7+1) M* 


(bio) 


(Bll) 


The ratio h m/R is plotted as a function of Mach number in figure 13 . 

Consider now a typical example of how the parameter P can be 
employed to guide asymmetric nozzle design. Assume a nozzle is designed 
with no reflected wavelets at a Mach number of 1.2; in this case, then, 

P equalB 0. The nozzle is required, however, to operate up to a Mach 
number of 2.0. It is clear from figure 13, the equation for P, a nd the 
previous discussion that, up to a Mach number of 1.8, P is less than 0 
and hence reflected compression wavelets are required near the exit. 

Above a Mach number of 1.8, P is greater than 0 and hence reflected 
expansion wavelets are required to give uniform flow at the exit. This 
change in the type of reflected wavelets required is difficult to obtain 
with walls of reasonable shape. Reflected compression wavelets result 
from an upper wall of relatively large curvature and reflected expansion 
wavelets result from an upper wall of relatively small curvature. 
Therefore, the design of a nozzle operating over a range of Mach num- 
bers spanning the peak in the h^/R curve is difficult. One nozzle 
was designed to give uniform exit flow at M = 1.2, 1.6, and 2.0; however, 
the second derivative of one portion of the upper wall varied rapidly 
over a short distance. This portion produced a shock wave at M = 1.8. 
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TABLE I.- COORDINATES OF NOZZLE W AT.TS DESIGNED FOR 
MACH NUMBERS OF 1.5 AND 2.6.* 


Xu 





% 

YU 




* L 

0.000 

0.0000 

9.200 

0.6860 

0.000 

0.0000 

6.500 

0.3629 

.200 

.0000 

9.400 

.7265 

.200 

.0000 

6.800 

.3976 

.400 

.0001 

9-600 

.7683 

.400 

.0001 

7.000 

.4345 

.600 

.0002 

9.800 

.6115 

.600 

.0003 

7.200 

.4737 

.600 

.0006 

10.000 

.8559 

.800 

.0007 

7.4oo 

.5153 

1.000 

.0011 

10.200 

.9016 

1.000 

.0013 

7.600 

.5592 

1.200 

.0020 

10.400 

.9486 

1.200 

.0022 

7.800 

.6056 

1.400 

.0032 

10.600 

•9969 

1.400 

.0035 

8.000 

.6544 

1.600 

.0047 

10.800 

1.0465 

1.600 

.0052 

8.200 

.7056 

1.600 

.0067 

11.000 

1.0975 

1.800 

.0074 

8.400 

.7593 

2.000 

.0091 

11.200 

1.1497 

2.000 

.0102 

8.600 

.8155 

2.200 

.0120 

11.400 

1.2033 

2.200 

.0134 

8.800 

.8743 

2.400 

.0154 

11.600 

1.2582 

2.400 

.0173 

9.000 

.9356 

2.600 

.0193 

11.800 

1.3146 

2.600 

.0219 

9.200 

.9995 

2.800 

.0238 

12.000 

1.3723 

2.800 

.0271 

9.400 

1.0660 

3.000 

.0290 

12.200 

1.4316 

3.000 

.0332 

9.600 

1.1352 

3.200 

.0348 

12.400 

1.4928 

3.200 

.0401 

9.800 

1.2071 

3.400 

.0412 

12.600 

1.5562 

3.400 

.0479 

10.000 

1.2818 

3.600 

.0484 

12.800 

1.6221 

3.600 

.0568 

10.200 

1.3592 

3.800 

.0563 

13.000 

I.6906 

3.800 

.0 668 

10.400 

1.4395 

4.000 

.0651 

13.200 

1.7618 

4.000 

.0780 

10.600 

1.5226 

4.200 

.0747 

13.400 

1.8357 

4.200 

.0905 

10.800 

1.6087 

4.400 

.0853 

13.600 

1.9123 

4.400 

.1044 

11.000 

1.6977 

4.600 

.0968 

13.800 

1.9914 

4.600 

.1197 

11.200 

1.7896 

4.800 

.1094 

14.000 

2.0730 

4.800 

.1365 

11.400 

1.8845 

5.000 

.1231 

14.200 

2.1568 

5.000 

.1548 

11.600 

1.9826 

5.200 

.1379 

14.400 

2.2423 

5.200 

.1748 

II.8OO 

2.0837 

5.400 

.1538 

14.600 

2.3295 

5.400 

.1963 

12.000 

2.1880 

5.600 

.1708 

14.800 

2.4183 

5.600 

.2196 

12.200 

2.2955 

5.600 

.1890 

15.000 

2.5087 

5.800 

.2445 

12.400 

2.4062 

6.000 

.2084 

15.200 

2.6006 

6.000 

• 2713 

12.600 

2.5201 

6.200 

.2289 

15.400 

2.6940 

6.200 

.2999 

12.800 

2.6374 

6.400 

.2507 

15.600 

2.7891 

6.400 

.3304 

13.000 

2.7581 

6.600 

.2737 

15.800 

2.8859 





6.800 

.2979 

16.000 

2.9642 





7.000 

.3233 

16.200 

3.0843 





7.200 

.3499 

16.400 

3.1860 





7.400 

.3778 

16.600 

3.2894 





7.600 

.4070 

16.800 

3.3945 

For M = 1.5, b = 2.236 | 

7.800 

.4345 

17.000 

3.5014 





8.000 

.4692 

17.200 

3.6100 


M = 2.6, b - 4.800 

8.200 

.5022 

17.400 

3-7204 





8.400 

.5364 

17.600 

3.8325 





8.600 

.5719 

17.800 

3-9465 





8.800 

.6087 

18.000 

4.0622 





9.000 

.6467 

18.200 

4.1798 






^Coordinates do not include correction for boundary layer 
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(a) - Physical plane . 


Mach angle, a M ,deg 
50 30 



(b) - Hodograph plane. 


Figure 3 -Typical construction of characteristic net at lower design 

M a ch number 
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90 70 50 30 G^H 

Mach angle t a M ,d eg 

(b)— Hodograph plane. 


Figure 4 . — Typical construction of characteristic net in reverse of 
flow direction at higher design Mach number . 


Streamline 
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Mach angle , a u ,deg 



8 



S3 

m 

HA.' 

vSffl. 

1 
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/ 
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(b) — Hodograph plane . 


Figure 5 — Typical construction downstream 

at higher design Mach 



of minimum section 
number. 
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Figure 6. — Combined constructions at higher design Mach number 
showing completed characteristic net and streamline uniting 
the two portions of the upper wait . 



Figure 7.- Construction for high Mach number design procedure. 
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Horizontal distance from leading edge 
of test section, test section heights 


(a) Along the center line. 



Figure 8.- The variation of the stream static pressure 
coefficient axially in the 2— by 2— inch test section. 





Stream pressure coefficient 
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Horizontal distance from leading edge 
of test section, test section heights 


(b) 0. 75 inches above the center line. 



Figure 8 — Continued. 





Stream pressure coefficient, 
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Horizontal distance from leading edge 
of test section, test section heights 


(c) 0. 25 inches below the center line . 



Figure 8. — Concluded . 
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A 



Figure 10.- Downstream end of flow field in an 
asymmetric nozzle with uniform exit flow. 



Figure II Section of asymmetric nozzle made up of 
elements of two Prandt! - Meyer streamlines. 






